Parity effects in the scaling of block entanglement in gapless spin chains 
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We consider the Renyi a-entropies for Luttinger liquids (LL). For large block lengths £ these are 
known to grow like In I. We show that there are subleading terms that oscillate with frequency 2/cf 
(the Fermi wave number of the LL) and exhibit a universal power-law decay with I. The new critical 
exponent is equal to K/(2a), where K is the LL parameter. We present numerical results for the 
anisotropic XXZ model and the full analytic solution for the free fermion (XX) point. 
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Luttinger liquid (LL) theory describes the low-energy 
(large-distance) physics of gapless one-dimensional mod- 
els such as quantum spin chains and correlated electron 
models. It corresponds to a conformal field theory (CFT) 
with central charge c = 1 and is known to provide accu- 
rate predictions for universal properties of many physical 
systems. LL theory has been applied successfully to re- 
cent experiments on carbon nanotubes [1], spin chains 
[2], and cold atomic gases [3]. A much studied example 
of a lattice model that gives rise to a LL description at 
low energies is the spin- 1/2 Heisenberg XXZ chain 

H = -X>Mw + «Kf i + aoKh] ■ (i) 

3=1 

Here (Xj are Pauli matrices at site j and we have im- 
posed periodic boundary conditions. Recent years have 
witnessed a significant effort to quantify the degree of 
entanglement in many-body systems (see e.g. [1] for re- 
views). Among the various measures, the entanglement 
entropy (EE) has been by far the most studied. By par- 
titioning an extended quantum system into two subsys- 
tems, the EE is defined as the von Neumann entropy of 
the reduced density matrix pa of one of the subsystems. 
The leading contribution to the EE of a single, large block 
of length £ can be derived by general CFT methods [SHI] • 
The case of a subsystem consisting of multiple blocks re- 
quires a model dependent treatment, but the EE can still 
be obtained from CFT [5]. On the other hand, little is 
known with regard to corrections to the leading asymp- 
totic behaviour. In the following we consider the Renyi 
entropies 

S a = - i —liiTrp1, (2) 
1 — a 

which give the full spectrum of pa [S] and are fundamen- 
tal for understanding the scaling of algorithms based on 
matrix product states [TUHH] ■ We note that Si is the von 
Neumann entropy and Soc gives minus the logarithm of 
the maximum eigenvalue of the reduced density matrix 
(known as single copy entanglement [131 El)- According 



to CFT, in an infinite gapless one-dimensional model a 
block of length I has entropies [SHI] 

^ CFT W = ^(l + ^)ln^ + C ' Q , (3) 

where c is the central charge and c' a a non-universal 
constant. In a finite system of length L, the block 
length € in (JSj) should be replaced with the chord dis- 
tance D(£,L) = ^sin^. In many lattice models the 
asymptotic scaling is obscured by large oscillations pro- 
portional to (— iy. Some typical examples are shown in 
Fig. [I] where we plot S a (i,L) for a = 1,2, oo for the 
XXZ model at A = —1/2 as obtained by density matrix 
renormalization group (DMRG) computations. While Si 
is smooth, is seen to exhibit large oscillations. For 

a = oo in particular it is difficult to recognize the CFT 
scaling behaviour While such oscillations have been 
observed in several examples jTS) [T6] and can be seen to 
arise from strong antiferromagnetic correlations, a quan- 
titative understanding of these features was until now 
lacking. We show that these oscillations obey the univer- 
sal scaling law 
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FIG. 1: Parity effects in Renyi entropies in the XXZ model 
at A = -1/2. S a (£,L) for several L and a = l,2,oo. The 
straight lines indicate the asymptotic slopes (1 + a~ )/6. 



2 



where p a is a universal critical exponent equal to 2K/a. 
Here K the LL parameter, kp is the Fermi momentum, 
and f a is a non- universal constant. In a finite system, 
the block length £ in Q is replaced by the chord dis- 
tance, and f a is multiplied by a universal scaling func- 
tion F a (£/L), that in general depends on the parity of 
L. We note that in zero magnetic field (half- filling) we 
have fc_p = 7r/2 and the oscillating factor in Q reduces to 
(— 1) as observed. While we establish (4) for the partic- 
ular case of the Heisenberg XXZ chain (1 ), where the LL 
parameter is given by K = tt/(2 arccos A), we expect the 
scaling form to be universal, because it is related to the 
low-energy excitations of the model and is therefore en- 
coded in the continuum LL field theory description. Re- 
cent results for the entanglement entropy confirm these 
expectations [T7] . 

XX model. This case corresponds to A = in ([T]). 
The LL parameter and exponent in Q are K — 1 and 
p a = 2/a respectively. The computation of the Renyi en- 
tropies can be achieved by exploiting the Jordan- Wigner 
mapping to free fermions, which reduces the problem to 
the diagonalization of an £ x £ correlation matrix (see 
[H] for details). Jin and Korepin (JK) showed [H] that 
Renyi entropies can be obtained by the following contour 
integral encircling the segment [—1,1] of the real axis 

S a (£)^^<[ea(l) dh ^ { ^dl. (5) 
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the Fisher-Hartwig formula was used to determine the 
asymptotic scaling of the Renyi entropy with £, which 
agrees with the CFT formula Here we employ the 
generalized Fisher-Hartwig conjecture |20j in order to go 
beyond the results of [19] and determine the subleading 
corrections. The terms in the asymptotic expansion of 
the determinant relevant for calculation of the Renyi en- 
tropy can be cast in the form 



Di(l) 
Df K {l) 
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where Dj K is the result of [TS], 2Tri/3{x) = ln[(l+x)/(l- 
x)], and Lk = \2£ sin kp\. The calculation of the integral 
in ([5]) is now straightforward. One expands In D# in ^ 
in powers of Lk, determines the discontinuity across the 
cut [—1,1], changes the integration variable from A to 
—i/3(l) and finally obtains the leading behavior from the 
poles closest to the real axis (details will be reported 
elsewhere [21). The resulting asymptotic expression is 
valid at fixed a, in the limit InLk 3> a. The final result 
is given by Eq. Q with 
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FIG. 2: Corrections to scaling d a (£) = S a (£) - S% FT (£) for 
the XX model and four different values of a. Red circles and 
black squares correspond to even and odd £ respectively. 



We note that /1 = and therefore no oscillating cor- 
rections for the Von-Neumann entropy are predicted, in 
agreement with numerical observations. 

The requirement that In 3> a implies that the 
asymptotic behaviour is only reached for very large block 
lengths, e.g. at a = 10 we need Lk ^> 20000. In the 
pre-asymptotic regime there are several sources of cor- 
rections. First the integral is no longer dominated by the 
poles closest to the real axis, which leads to power law 
corrections of the form ^^ 2m / a ( w jth integer m), which 
oscillate as e ±l2k f e . Corrections with different oscilla- 
tory behaviour arise from the higher order terms in the 
expansion of lnD^(/) in powers of L^. The first correc- 
tion is proportional to e ±tik F l ; the next to e ±t6kFe etc. 
In zero magnetic field, where kp = tt/2, the leading term 
is proportional to (— ) e while the second does not oscil- 
late. Hence there is a subleading constant background 
in addition to the leading oscillatory behaviour. In the 
limit a — > 00 all these terms become of the same order, so 
that we need to resum the entire series that arises from 
expanding lnZ)^(Z) and then carrying out the A integral. 
In the zero magnetic field case we thus obtain 
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Here the constant b w 7.1 has been fixed by summing 
certain contributions to all orders in 1/ihiLk) and agrees 
well with numerical results 22] . 

Numerical results for the XX model can be obtained 
by diagonalizing the correlation matrix both infinite and 
finite systems. We first present the results for infinite 
systems. We consider only the model in zero magnetic 
field and plot the quantity 



da(£) = S a (£)-S^ L (£) 



(9) 
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FIG. 3: Corrections to scaling in the finite length XX model 
for odd L. F2(x) and Foo(x) (x = £/L) obtained from 69 
different systems with lengths in the range 17 < L < 4623 
exhibit perfect data collapse. 



where the value for the constant contribution cL in 
S% FT {£) is taken from 19 . According to Eq. Q for 
k F = tt/2, d a {i) ~ {-) e (2e)-P-f a . In Fig§we compare 
the absolute value of d a (£) for a = 2,5, 20, oo and block 
sizes £ up to 4000 sites to our asymptotic results Q , (|8]) . 
For a = 2 the curves for odd and even £ are practically 
indistinguishable (the line corresponding to the analyti- 
cal result is invisible under the data points). For a = 5, 
we still obtain power laws with the exponent p 5 = 2/5, 
but the curves are not as symmetrical as for a = 2 be- 
cause subleading corrections become visible. Increasing 
a further, the deviations of d a {l) for £ < 4000 from the 
asymptotic behaviour become quite pronounced. This 
is shown in Fig|2] for the case a = 20, where the lead- 
ing asymptotic result (straight line) is seen to be a poor 
approximation to d a (l) for even £. Including the sublead- 
ing corrections gives curves perfectly covered by data in 
Fig(2] The last panel in Fig(2] shows the result in the 
a = oo. The numerical results are seen to be in perfect 
agreement with Eq. ([8]). 

We now turn to finite systems. We numerically deter- 
mined the quantity (recall that D(£, L) = ^ sin ^) 

F a (£/L) = [S a (£,L) S^ FT {i,L)]f- l D(i,Lf' a , (10) 

(for a — oo we multiply by hibD(£, L)) for a variety of 
values of a and system sizes ranging from L — 17 to 
L = 4623. We observe that there is data collapse for 
any L on two scaling functions for £ odd and even re- 
spectively. Results for the cases a = 2, oo and odd L are 
shown in Figj3j The quality of the collapse is impres- 
sive considering that there are no adjustable parameters 
and that the plots contain millions of points ranging over 
three orders of magnitude in both £ and L. For a=2we 
observe that F%{x) — ±cos7rx (these are shown as con- 
tinuous curves in Figj3]). We currently have no analytical 
derivation of this scaling function. For other values of a 
we obtain similar data collapse, but the quality decreases 
with increasing a, indicating the presence of other cor- 
rections. For even L we obtain different scaling functions 
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FIG. 4: Left: The additive constant in the Renyi entropies c' a 
(bottom) and the multiplicative constant c a for the moments 
Trp5. (t°P in log-scale) as function of a at fixed A. The drawn 
lines are the exact values for A = 0. Right: c' a as function of 
A at fixed a. 



- F a (x) then is almost constant (see below) . 

XXZ model and DMRG. To characterize the XXZ 
model in the gapless phase with — 1 < A < 1 we per- 
formed extensive DMRG calculations at finite L. We 
used the finite-volume algorithm keeping x — 800 states 
in the decimation procedure. This rather large value of 
X is required to obtain a precise determination of the full 
spectrum of the reduced density matrix in the case of 
periodic boundary conditions. The data we have used 
in our analysis can be considered as numerically exact. 
Hence the main limitation as compared to the XX case is 
the relatively small value of L accessible by DMRG (we 
considered 21 < L < 81 for odd L and 20 < L < 80 for 
even L). Another complication stems from the fact that 
the value of the constant contribution c' a to the Renyi en- 
tropy is not known analytically for A ^ 0[53]. We obtain 
it by fitting our numerical data. The results are shown 
in Fig. [4] The data for A = is in good agreement 
with the exact results of [19] (full line), establishing the 
correctness of our fitting procedure and the reliability of 
DMRG. The multiplicative constant c a = e( 1 ~ a ^ c <* in the 
moments of the density matrix Trp^ = Ca^/ 6 *- 0-1 /") in 
Fig. [4] shows an exponential decay with a, except for a 
very close to 1. Hence c a can essentially be absorbed into 
a rescaled block length £ as was pointed out in [9l [23] ■ 

Having obtained the constant contribution d a one can 
determine the universal scaling functions F a (x). We 
present results for both even and odd L and a number 
of representative values of a and A in Figj5] We find 
that the data collapse is very good for all cases. This 
is remarkable given the limited system sizes accessible 
by DMRG. We note that, as expected, the data collapse 
becomes poor in the vicinity of the two isotropic points 
A = ±1. At A = — 1 there is a marginal operator (see e.g. 
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[T5] ) that gives rise to well known logarithmic corrections 
to scaling for correlation functions. In the ferromagnetic 
limit A — >• 1 the model loses conformal invariance (the 
dispersion relation becomes quadratic) and is no longer 
described by a LL. Hence none of the results presented 
here is expected to hold. 

The critical Ising chain has no strong antiferromag- 
netic correlations and we therefore expect the corrections 
to scaling to be non-oscillatory. It is easy to see that this 
is indeed the case. Igloi and Juhasz [21] showed that the 
EEs of the XY chain can be expressed in terms of the 
EEs of two Ising chains. At the quantum critical point 
this relation reads S* X (2£,2L) = 2S l a {l,L), where S\ 
refers to the critical Ising chain. This implies that the 
Renyi entropies in the Ising chain are just one half of the 
corresponding entropies in a XX chain of twice the block 
length and twice the system size. As both I and L are 
even, our results for the XX model imply that the correc- 
tions to scaling are non-oscillatory and decay as l~ 2 l a . 
This agrees with numerical computations. 

Conclusions. In this letter we considered the Renyi 
entropies for the critical spin- 1/2 Heisenberg XXZ chain 
with periodic boundary conditions. By a combination 
of analytic and numerical techniques, we computed os- 
cillating corrections to scaling which are expected to be 
universal. These are parametrized in terms of the Lut- 
tinger parameter K and the Fermi momentum kp. We 
argued that our results hold generally for Luttinger liq- 
uids. In the case of open boundary conditions, a similar 
relationship is expected to hold, but with the replace- 
ment K — > K/2 (see also [H]). It would be interesting 
to prove this at least in the special case of free fermions. 
Finally we would like to comment on our results in light 
of a recent proposal [5], that one way of distinguishing 
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FIG. 5: Universal scaling function f a F a (x) (where x = t/L) 
for the XXZ chain for several values of a and A. Left: Several 
a and A for odd L. Right: even L. In the latter case F a (x) 
is practically independent on x. 

between different theories with the same central charge 



is to consider the entanglement of multiple blocks. Our 
results establish that it is sufficient to consider a single 
block once one takes into account the universal sublead- 
ing oscillatory corrections. 
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